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1. Introduction 

Despite numerous efforts conducted in the past, the understanding of macroscopic 
systems out of equilibrium is still far from being put in a systematic theory. Unlike 
the relaxation around the equilibrium which has become a standard textbook material, 
less can be said for systems that are maintained far from equilibrium. One possible 
route to fill this gap that proved useful in the past (e.g. in developing general theory 
of critical phenomena) is to study particular microscopic models. Usually, one starts 
with a minimal model that admits analytical treatment (exact or approximate), later 
adding more details to make it more realistic. Unfortunately, the lack of detailed 
balance - a defining property of systems maintained far from equilibrium - means that 
we are generally deprived even of the knowledge of its steady state, not to mention 
the relaxation mechanism towards it. Exactly solved models far from equilibrium are 
thus rare, but have an important role in nonequilibrium statistical mechanics. 

One such model is the asymmetric simple exclusion process (ASEP), a minimal 
model of transport of (classical) particles driven by an external field and interacting 
only through the exclusion principle that prevents them from coming too close to each 
other. Although simplified, this interaction describes several real situations on various 
length scales, ranging from mobile ions in superionic conductors (T] to self-propelled 
particles in mesoscopic (ribosomes on a mRNA [3J|3]) and macroscopic (cars [I]) 
systems. From a theoretical viewpoint, the ASEP has become a paradigmatic model 
of boundary-induced phase transitions that are, unlike their equilibrium counterpart 
(with short-range interactions), present even in one dimension. Originally proposed 
to model translation of mRNA more than five decades ago (2j|3] , the exact solution in 
1993 [5,6 sparked a great interest that led to several important results relevant to the 
general theory of nonequilibrium steady states. Using the exact solution of the ASEP, 
Derrida, Lebowitz and Speer [7] derived a non-equilibrium analogue of the free energy 
that describes coarse-grained fluctuations around the steady state, thus generalizing 
older Onsager-Machlup theory j8] valid around equilibrium. This later helped Bertini 
et al [9] to build the macroscopic theory of fluctuations for driven diffusive systems, 
of which the ASEP is just one example. 

A lot of work has been devoted to improve the ASEP to better fit particular 
phenomena, most of which have their origins in vehicular traffic or biology (for a 
comprehensive and recent review see |10j). The exact solution has been successfully 
extended to particle-dependent hopping rates |lHfl"4] and multispecies systems [15lfT6] , 
both of importance for traffic phenomena. In biology, some of the generalizations 
include particles occupying more than one site [TTJUH], position-dependent hopping 
rates (termed site- wise disorder) [19] . desorption and adsorption of particles in the 
bulk (i.e. Langmuir kinetics) |20] , particles with internal states |21) , extension to more 
than one lane [22] . dynamically extending lattice [23] , etc. Unfortunately, common to 
most of these generalizations is non-applicability of the exact methods that were used 
to solve the original ASEP, namely the matrix-product ansatz (see e.g. [24]). Instead, 
most approaches utilize various mean-field approximations that neglect correlations 
in a fashion similar to truncation in the Bogoliubov-Born-Green-Kirkwood-Yvon 
(BBGKY) hierarchy. While they give a satisfactory account for most of the 
aforementioned phenomena, it is less clear how to improve them in some controlled 
fashion. In some cases such as site-wise disorder, the mean-field approximation itself 
becomes hard to treat analytically as the number of inhomogeneities increases. 

In this article we study totally asymmetric simple exclusion process (TASEP) 
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- a version of the ASEP with unidirectional hopping - in the presence of site-wise 
disorder and show that some of the information about its exact steady state can still 
be retained. Our starting point is the known fact, reviewed in |Appendix A| that 
the non-normalized steady-state weights which solve the master equation can always 
be written as multivariate polynomials in all the hopping rates or alternatively, as 
univariate polynomials in one of the hopping rates with polynomial coefficients that 
depend on all other hopping rates |25) . As our main result, we show how to compute 
the first few lowest polynomial coefficients exactly, which gives a good approximation 
of the exact steady-state weights if the hopping rate that we expand in is much smaller 
than the other rates. Since in that case the particle current is small, we term this 
regime the low-current regime. Our method of computing non-normalized steady- 
state weights in the low-current regime is not restricted to site-wise disorder but can 
be readily applied to many other generalizations of the TASEP and other driven 
diffusive system with unidirectional hopping. 

The paper is organized as follows. In section [5] we define the TASEP with site- 
dependent hopping rates and review some of the related exact results. A general 
approach for calculating steady state in the low-current regime is devised in section 
[3] for arbitrary disorder distributions. Section |4] is devoted to few particular cases of 
binary disorder where all the slow sites share the same hopping rate r < 1. Our main 
interest there is to calculate steady-state average of the particle current, i.e. its Taylor 
expansion in one of the hopping rates which is considered to be small. Some of the 
exact results that we obtain have already been known, but here they are derived for the 
first time rather then being guessed from the exact solution of small systems |17 [ l26 j . 
Further applications, mainly devoted to the most general case of non-binary disorder, 
are discussed in [5] with some interesting new implications for the protein synthesis. 

2. Model 

We consider the totally asymmetric simple exclusion process (TASEP) on an one- 
dimensional lattice of L sites, where each site i — 1,...,L is either occupied by 
a particle (rj = 1) or empty (rj = 0). Particles each move forward stochastically 
at rate Pi (which is modelled by the random-sequential update), provided the site 
i + 1 in front is empty (exclusion principle). In this paper we consider mainly binary 
disorder [371[5S], which means that Pi is either r < 1 (slow sites) or 1 (regular sites). 
(Non-binary disorder with pi 6 R to which our approach applies as well is left for 
section[5]) Boundary conditions can be either periodic (tl+i = t±) or open, the latter 
meaning that new particles are injected at site site i — 1 at rate a provided it is empty, 
and are removed from the site i = L at rate (3. An illustration of the process (not 
drawing boundary rules explicitly) is presented in figure [T] 

r l r 




Figure 1. Schematic picture of the TASEP with binary disorder. Bonds where 
particles jump at rate r < 1 are represented by (permeable) walls to emphasize 
that letting r = means no particle is allowed to move through the wall. 

The steady state is described completely by the probability P{C) to find 
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the system in a given configuration C — {ti, . . . , tj,}, which solves the following 
continuous-time master equation, written here in a compact form, 



(periodic) = J^Pijn - r i+ x)P(. . ■ , Tj-i, 1, 0, n +2 , ■ ■ ■) (la) 



(open) = -a(l - 2 n )P(0, . . .) + (3(1 - 2r L )P(. ..,!) + 



L-l 



+ ^2pi{n - n + i)P(. . . ,^-1,1,0,^+2, . . .)• (lo) 

i=l 

In | Appendix A| we show that the steady-state probability P{C) is proportional to 
the determinant of a matrix whose matrix elements are linear combinations of hopping 
rates. That allows us to write P(C) in the following form, 

where /(C) is a multivariate polynomial in all hopping rates, 

(periodic) f{C) = Y,fk{C)r k (3 a) 

k 

(open) f(C) = J2kjAC)a l Pr k . (36) 

The particular form of ([2]) with /(C) given by (f3~aj) or (|3~6jl means that the ensemble 
average of any physical observable g(C) (e.g. local density Ti at site i, or current 
PiTi(l — i~j+i)) is a rational function with respect to any of the hopping rates (say r) 



(a(O) - V a(C)P(C) - 52c 9(C) f{C) _ E fc ^r k 
(g(C))-^g(C)P(Q- Ecf{c) - Ekhrk ; 



(4) 



c 

where coefficients a k and b k are given by 



a k = Y,g(C)f k (C), b k = J2fk(C). 
c c 



By expanding in Taylor series around r = we arrive at 

oo 

( 5 (C))=^c fc r fc , (5) 



fe=0 



where 



— -2^c fe _„— , fc = o, 



6o ^ 6o 



-Vc fc _„-^, k>M + l. 

n—l 

Here A4 is the maximal degree of all /(C)'s as polynomials in r, and can be calculated 
by invoking Schnakenberg's network theory [29) . The value of M. is not important 
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in our approach, as calculating coefficients beyond first few terms in ((5]) becomes 
extremely difficult. 

The idea of expanding (g(C)) in small r originates from the work of Janowsky 
and Lebowitz |26j . They obtained the first few coefficients in the expansion of current 
Ji(r) around r = in the TASEP with a single slow site both in the periodic and in 
the open boundaries case with r <C a — (3 

3 2 19 3 21535 4 77729356627 5 rw 6 . , . 

j L (r)=r r 2 H r J r H r 5 - 0(r b ). (6) 

J w 2 16 27648 146767085568 v 7 w 

The expansion (|6|) was obtained by solving the full master equation for small systems 
with L < 8 and noticing that as L increases, low-order terms become independent of 
L. A similar approach, coined finite segment mean-field theory (FSMFT), has been 
devised by Chou and Lakatos in [TH] for the open boundaries case with few slow sites. 
When slow sites are confined to a small segment, the idea is to solve the master 
equation exactly for that small segment and then calculate the current j(er_,<7 + ) as 
a function of the average densities cr_ and a+ at its ends. Assuming that the density 
profiles are flat outside the segment, <r_ and er + can be then calculated numerically 
by solving <r_(l — cr_) = j(<7_,<7 + ) = <r + (l — a + ). For r <C a, (3 > 1/2, Chou and 
Lakatos were able to deduce the analytical expression for the current up to 0(r 2 ) for 
two particular configurations of disorder, namely 

jL{r) = dTT r + 0(r2) ' (7) 

for two slow sites placed d sites apart, and 

3L(r) = ±±±r + 0(r 2 ), (8) 

for a bottleneck of I slow sites. As d -> oo in the current for r <C a, (3 > 1/2 
approaches Jl(V) k r + 0(r 2 ) as in the TASEP with one slow site. In ©, the current is 
clearly dominated by the capacity of the bottleneck and approaches r/4 when I — > oo. 
The downside of FSMFT is that the present computing power restricts the size of the 
segment that is treated exactly to ps 20 sites, thus excluding more complex disorder 
configurations. Also, FSMFT is basically a brute force attack on the master equation 
and tells us little about where the coefficients e.g. in <JHJ)-((H1) come from. In the 
next section we will present a general approach for computing low-order terms in the 
expansion f(C) — /o(C) + 7/1 (C) + 0(7 2 ), where 7 stands for any of the model's 
hopping rates such that the current ^(7) — > as 7 — s> (7 = r in the periodic 
boundaries case or 7 € {r, a, ft} in the open boundaries case). 

3. Main idea and general results 

We start by writing /(C) as a polynomial in one of the model's hopping rates 7, 

Z(C)=$>(C) 7 fc , (9) 



where we assume fk(C) to be dependent on C and all other hopping rates 7^ 7. 
Inserting © in the stationary master equation (|lal) or (fT"6t and collecting all the 
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terms of 0(7 fe ) gives recursion relations 

= - 5 7iWi c^c))W(C' C)/ (C") 

c ' (10) 
- 5> " <W(c^c)MC -> C")/o(C) 
c 

= ^(1 - S l!W(c ^c))W(C' -> £7)/ k (C') + J] « 7 ,^(C'-,o))/fc-i(C) 

c ' c ' (11) 

- 5> - Sy,w(c^c>))W(C -y C")/ fc (C) - XXvrCc-^oJ/fc-iCC). fc > 0. 
c c 

It is useful to picture our system as a directed graph made of vertices (configurations) 
and directed edges (transitions between configurations) weighted by the hopping rates. 
Edges that are weighted by 7 are called slow, and all the others are called regular. A 
special role here is played by configurations that have all their outgoing edges slow. We 
will call such configurations blocked (with respect to 7) because if the system gets into 
one of these configurations and 7 — > 0, it will stay there forever. If B is a non-empty 
set of all such configurations then terms /o(C) are clearly absent from (|T0|) for any 
C £ B. For example, the TASEP with a slow site and periodic boundary conditions 
has only one blocked configuration with respect to r, the one in which all particles are 
behind the slow site. In the open boundaries case there are two additional blocked 
configurations, one with respect to a (an empty lattice) and one with respect to /? 
(full lattice). 

The fact that some terms are missing when collecting zeroth-order terms may 
seem confusing, because we have to go to the next order in 7 to calculate fo{C), but 
we need /o(C) to calculate the next order terms. In the next section we show how to 
eliminate all first-order terms, thus ending up with a closed set equations for fo(C), 
CeB. 



3.1. Zeroth-order terms 

We start by outlining the general procedure for calculating zeroth-order terms and 
then show how it works on an explicit example. Since /o(C) is absent from (|10[) for 
any C £ B, we can assume that /o(C) 7^ for C £ B and then solve (fit))) by setting 
all the remaining zeroth-order terms to 0, 

/o(C)=0, CiB. (12) 

Thus the existence of blocked configurations reduces the calculation of fo(C) to a 
smaller set of configurations B. 

To get a closed set of equations for the remaining zeroth-order terms we have to 
inspect terms that are linear in 7, 

=5^(1 - 5^w ( c^c))W{C -> C)h{C) + ^ S y ,w{c^c))fo(C) 
c c 

c c 



(13) 



It will prove useful to write (|13[) in the following form, 

h{C)= 53 A(C,C")/o(C"), C£B, (14) 

C'EB 
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where A(C, C) is some unknown matrix. Inserting (TLI)) in (fT3)) . the equation for fo(C) 
now reads 



Ay(C)f (C)= ^ k{C,C")MC% CeB (15) 
C"eB 

where A 7 (C) and k(C, C") are given by 

£ 7 ,VK(C-K7')> ( 16 ) 



C" 



k(C, C") = V(l - <y 7 , W (c^o)A(C", C") + 5^ w{c „^ C y (17) 



E 

c 

To find / (C), we have to find A(C,C*'), calculate k(C",C") and then solve (fT5|). 
Luckily, k(C , C") can be found without using the expression (|T7|) . The algorithm 
that we give below is essentially the same as the one we'll use later for constructing 
A(C, C). 

To start with, let's call a path V in configuration space any sequence of 
configurations Ct,C%,. . . ,C„ such that none of W(C\ — > C2), W(C2 — >• C3), . . ., 
W(C n -i — > C n ) are zero. A regular path is a path in which none of the edges is slow. 
With all these preliminaries, we rewrite the equation (|13l) for fo(C), C G B, 

Ay(C)f (C) = W(C -> C)(l - £ 7 

C ' , (18) 

+ E ^r,w(Cf'->cr)/o(C")i CeB. 
c 

Configurations C" on the r.h.s. are obtained by moving one of the particles backwards. 
Any move to a configuration C ^ B across a slow edge should be discarded since 
fo(C) = for C £ B and the second term on the r.h.s. vanishes in that case. Now 
let's focus on C for which W(C — >• C) 7^ 7. Since we started from CeB and 
moved one particle backwards across a regular edge, C has the following properties: 
(a) C £ B, (b) /o(C") = and (c) -4 (C") = 1. The last one is simply because C 
is just one jump from C and therefore cannot have any other regular outgoing edges 
except the one pointing towards C. Now, let's write the equation for /i(C"), 



MC) ■fi(C') + A r (C')f (C) =J2 S -y,w(C"^cnfo(C") 



=0 



+ V W(C" -> C")(l - * 7 ,w(o»-.oo)/i(C")- 



(19) 



C" 



The l.h.s. of (fT9"|) is exactly what we get on the r.h.s. of (fTE)) by moving one 
particle across a regular edge from CeB. We can then insert (| in (fTE)) and thus 
eliminate fi(C'). The idea is to repeat this process of moving particles backwards and 
substituting f\{C") from 

Mc")h{c") + Mc")h{c^ = J2 Kw{c»>^c>'MC") 

(20) 




+ £ -> C")(l - 5 1>w{cl ,^ cll) )h{C" 

C" 
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for any C" B that we reach along. Mathematically speaking, we must exhaust all 
backward paths originating from C that have one slow edge at the end and all other 
edges regular. Moreover, since /o(C') = for all C ^ B, we should consider only 
paths that end in configurations belonging to B. 

Now, consider any C" ^ B that is reached by moving particles backwards from 
C without crossing a slow edge. Then if we start at C" and move particles forward 
without crossing a slow edge, we must end at C. This ensures that by moving particles 
backwards in order to eliminate any fi(C") along the way, /i(C") will appear exactly 
Ao(C") times. We can then complete the l.h.s. of O (since Ay(C")f Q (C") = 0) 
and substitute (|20|) in (ITHl) . In the end, when all paths have been exhausted and all 
first-order terms eliminated, the final result is 

A y (C)/ (C)= ^ C '^ CeB ( 21 ) 

c"eSo(C) 

where So(C) is the set of all configurations C such that (a) C can be reached from 
C by moving particles backwards and crossing a slow edge in the last move only and 
(b) 'fo(C') ^ 0. (From this definition, S Q (C) C 5 for any C £ B.) Going back to (JT7J) 
we have therefore proved that 




Figure 2. A part of the graph representing the TASEP on a ring of L = 5 
sites with N = 2 particles and slow sites placed at i = 2 and i = 5. Blocked 
configurations are designated with ellipses and slow edges with double arrows. 

To illustrate this, let's consider the TASEP on a ring with L = 5 sites, N = 2 
particles and two slow sites placed at i = 3 and i = 5. There are 5!/(2!3!) = 30 
configurations in total of which | B | = 3 are blocked, B = {o • • j o o | , o o • | o • | , o o o | • • | } 
(as in figure [TJ sites in front of the walls are considered slow). A portion of the graph 
containing configuration C = oo»|o»|is presented in figure [5] The equation for 
/ (o o«| o »|) reads 

2/ (o o • o »|) = /i(o • o| o »|) + fx(o o • • o|) = 

= LM° • °l ° •!) + /o(° • °| ° •!)] + /i(° ° H • °|), ( 22 ) 

s > 

^0 
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where in the second line we added / (o • o| o «|) = to complete the master equation 
for /x(o • o| o •!), 

f x (o • o| o »|) + / (o • o| o •!) = f t (m o o| o »|) + f^o • o| • o|). (23) 
Master equation for the remaining term /i(o o • • o|) reads 

/ 1 (oo«|.o|) = / 1 (o.o|«o|). (24) 
Substituting (|2"3"|) and (|24l) in (|2"12j) . the equation for /o(o o • o »|) now reads 

2/ (o o . o .|) = o o| o .|) + 2/i(o • o| • o|). (25) 
Equations for the terms on the r.h.s. of (f2"5"j) are 

/!(• o o| O «|) = o o| • o|), 

2/i(o • o| • o|) = o o| • o|) + / (o • • o o|), 
which upon substitution in (|25p becomes 

2/ (o o . O .|) = 2A(. O o| . o|) + / (o . . O o|). (26) 
Finally, inserting 

2/i(» o o| • o|) = / (» o • o o|) +/ (o o o| • •!), 

v , ' 

=0 

in (1261) gives the final equation for / (o o »| o »|) 

2/ (o o • o »|) = / (o • • o o|) + / (o o o| • »|). (27) 

This daunting task that we have just performed in fact has a remarkably simple 
interpretation. To see it, let's write equations for the remaining blocked configurations, 

/ (o..|oo|)=/ (oo.|o.|) ! (28) 
/o(ooo|..|) = / (oo.|o.|). (29) 

The process described by (|2"7|) -(f2"9 |) can be interpreted as follows: any particle that 
jumps from a slow site immediately joins the queue in front, while at the same time 
the whole queue that the particle has just left moves one step forward. Intuitively, this 
is easy to understand: the limit r — > creates a huge separation of time scales, so that 
any particle that jumps from a slow sites is likely to join the queue in front before 
any other jump from a slow sites happens. Particles are thus exchanged between 
compartments separated by walls, each compartment having a finite capacity which 
is the number of sites between two neighbouring walls. This number can be greater 
than 1, as in our example above, and so the process that TASEP maps to in the limit 
r — > is a generalization of the exclusion process called the partial exclusion process 
(not to be confused with partially asymmetric exclusion process). Interestingly, the 
partial exclusion process has been introduced long time ago by Schiitz [30], but has 
been rarely studied since [31"1l32]. 

For zeroth-order terms, our result can be summarized as follows: consider the 
TASEP with periodic boundary conditions on a lattice of L sites, N particles and D 
slow sites placed at sites kx, k-z, ■ ■ ■ ,kn — L (due to periodic boundary conditions, we 
have a freedom to place the last slow site at the end). Then /o(C) for C E B solves 
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always full always empty 



Figure 3. Mapping the TASEP with binary disorder in the limit r — > (left 
column) to a partial exclusion process (right column) in (a) periodic boundaries 
case and (b) open boundaries case. In both cases particles move to the right. 



the master equation of the totally asymmetric partial exclusion process on a lattice 
of D sites, where each site i = 1, . . . , D can hold at most ki — ki_\ particles, where 
fco = 1 (see figure [3h). 

One striking thing about this result is that k(C,C) and therefore any fo{C) 
are completely independent of the other hopping rates (if they exist) and instead 
depend only on the "connectivity" of blocked configurations. We will go back to this 
observation more in section [5] where we discuss two or more different slow hopping 
rates. For now, let's just focus on what this means for the binary disorder in the open 
boundaries case. In the open boundaries case, all blocked configurations (with respect 
to r) have the first compartment occupied by particles and the last one empty. If we 
have a lattice of L sites and D slow sites placed at hi, ■ . ■ , fcr>, then /o(C) for C e B 
solves master equation of a totally asymmetric partial exclusion process on a lattice 
of D — 1 sites, where each site i = 1,D — 1 holds at most ki+\ — ki particles. At 
the boundaries, particles jump into the first compartment at rate 1 if it's not full and 
leave the lattice from the last compartment at rate 1 (see figure [3]d). Thus the rates 
at which particles are exchanged with reservoirs are always maximal, i.e. equal to 1. 

For a small number of blocked configurations we may hope to solve (p?Tj) 
analytically (as in section 0} or numerically (even for L large) . As the number of 
blocked configurations increases, since the exact solution for the steady state of the 
totally asymmetric partial exclusion process is not known, we may end up with a 
problem no less harder that the one we started with. In the next section we give a 
general recipe for how to calculate first-order terms /i(C) if we can somehow solve 



3.2. First- order terms 

To calculate fi(C), C B, we have to determine X(C,C) in (Tl4l) . The equation for 
/i(C) for any C <£. B reads 

f 1 ^) = , rp^. y~] <5 7 ,w(C"^c)/o(C") 

+ — J^c) — K l ~ 5 i,w(c'^c))h{C ), 
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(32) 



where in the first sum we have used the fact that fo(C") = for all C" ^ B. Leaving all 
zeroth-order terms intact, we insert expressions like (|30[) recursively for the remaining 
first-order terms, until we are left with zeroth-order terms only. Again, we end up with 
as many zeroth-order terms as there are blocked configurations that can be reached 
from C by moving particles backwards but crossing a slow edge in the last move only. 
Recalling that the set of all such configurations is S'o(C), let's index a backward path 
from C to any C G Sq(C) with Vc,C> = C,C\,C%, . . . ,C n ,C (n can, of course, vary 
from path to path). For any given C, A(C, C) is then given by 

\<r r>) - V w{Cl c) w{ ° 2 ± c ^ w(c n - 2 -> Cn-j) 1 

^ MC) MOi) "' MCn-i) MCnV [ ' 

1 C,C 

for C G So (CO and is for C" $ S Q (C). 

To calculate fi(C) for C e B, we use the same recipe as for /o(C). Starting from 
C £ B, we look for paths Vc,C from C to any C' such that (a) C can be reached from 
C by moving particles backwards and crossing a slow edge in the last move only and 
(b) fx{C) 7^ 0. Let's call S\(C) the set of all such configurations and V\{C) the set 
of all configurations that are visited in going from C to all C £ Si(C) (not including 
C and C). Starting from the equation for /i(C), C £ B, 

Ay{C)h{C) = ]T W(C" -> C)(l - * 7)W (o^o)/a(C") 

C" 

+ E<Wc»^c)A(C"), Ce£, 

C" 

the idea is to eliminate W(C" — > C)/ 2 (C") by noting that 

MC^MC") + Ay{C")h(C") =Y,8 1<w{c „,^ c „ ) h(C'") 

+ J2W(C" -+ C")(l - <W(c»^c''))/2(0, 

C" 

This time, however, fi(C") on the left is not necessarily zero. To eliminate /2(C") 
in j32j, we add ^ 7 (C")/i(C") to both sides of ([32]) and substitute A(C")/2(C") + 
A 7 (C")fi(C") with (j3"3")l . By repeating the process of moving particles backwards and 
eliminating any Ao(C")f2(C) by adding and subtracting A^(C")fx(C"), we finally get 
a closed system of equations 

A 1 {C)h{C)= J2 E MC")h{C"), CEB. (34) 

Since S\(C) also contains blocked configurations, we can rewrite (|3"4")l as 

^ 7 (C)A(C)- ^ f 1 (C') = h 1 (C), CeB, (35) 
ceSi(C)ns 

where h\(C) is given by 

MC) = ^ C '^ E ^7(0/1(^0, (36) 

CeSi(C) C'eVi(C) 

C'gB 

If it weren't for the hi(C), would be just the same as in (|2"Tj) . Because of the 
non-zero terms on the right, (1331) no longer describes a stochastic process as in (f2"T|) . 
It may also be much difficult to solve (j34l) than (|2"Tj) if _B is large. 



(33) 
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3.3. Higher-order terms 

In principle, the same procedure can be applied to higher order terms. Let's denote 
with Sk(C) the set of all configurations C such that (a) C can be reached from C 
by moving particles backwards and crossing a slow edge in the last move only (b) 
/ fe (C) ^ 0. Then for C £ B we have 

f k {C)= J2 HC,C')f k - l (C')-^Qf k - 1 (C) ) C£B(37) 

where A(C,C") = for C £ S k -x(C), and is given by for C e S k -i(C). The 
additional term on r.h.s. of (|37|) was not present for k — 1 only due to the fact that 
fo(C) = for C i B. 

Now, for C £ B, let's denote with Vk(C) the set of all configurations that are 
visited in going from C to all C G S k (C). To find f k (C), C 6 B, we have to solve the 
following system of equations, 

A 1 (C)f k (C) = Y, E MC")h(C")> CeB. (38) 

Ces h (C) C"ev k (c) 

As for fc = 1, we can rewrite (|38p as 

Ar(C)MC)- fk(C) = h k (C), CeB, (39) 

C'eS fc (C)ns 

where h k (C) is given by 

fck(C9= E E AA(C'), CeB. (40) 

Ces k (c) C'GVfc(C) 

For the reasons evident in the following section, going beyond linear order becomes 
highly non-trivial for larger systems. In the rest of this paper we consider therefore 
only zeroth- and first-order terms in some simple disorder configurations. Insight that 
this approach gives us for general configurations of disorder is discussed in section [SJ 



4. Examples 

4-.1. TASEP with a single slow site 

Periodic boundary conditions. The TASEP on a ring with a slow site placed at site 
L has only one blocked configuration Cp with respect to hopping rate r < 1, the one 
in which all particles are immediately behind the slow site. It follows then that the 
equation (|38|) is trivially solved for any k > 0, i.e. f k {Cp) = const. = 1 for any k > 0. 
Using the fact that fo(C) = for any C £ B = {Cp}, we can write /o(C) using the 
delta Kronecker function, fo(C) — 5c.c P - To calculate the small-r expansion of the 
current iz,(r), we can choose jh = ( rT L{C)[l — ri(C)]), so that 

JiM^r-^ + OtA (41) 
bo b Q 

where a k and b k are given by 



a k =Y^(C)[l-n(C)]f k (C), b k = Yfk(C), k>0. 
c c 
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From here it follows easily that a = b = 1 and therefore jz,(r) = r + 0(r 2 ). To 
calculate the second-order term in r, we have to find fi(C), i.e. A(C, Cp) because of 

h(C) = HC, C")/o(C") = X! A ( c > CO^'.Cp = A(C, Cp). 
c c 

The construction of A(C, C"), as explained in section 13.21 tells us to look for 
configurations C such that the blocked state Cp is reached by moving particles 
backwards from C, provided the slow edge is crossed only by entering Cp. It is 
easy to see that any such C must have a particle displaced from the queue (figure 0]). 
Thus the configurations giving non-zero first-order terms are particle-hole excitations 
of the blocked state Cp . 




Figure 4. Configurations that have non-zero first-order term /i(C) ^ have 
either (a) one particle outside the queue and one hole inside the queue or (b) one 
particle taken out of the queue and added to its end. 



To calculate a\ and b±, it proves useful to define a set of configurations J7l, 

j L = {C\T L (C) = i, n (C) = o}. 

In other words, Jl is simply the set of all configurations such that rr^C^l — T^C)] ^ 
0. Using this definition, second-order term in Pdj) can be rewritten as 

c 2 = -^i = - £ X(C,Cp). 

This expression slightly simplifies the calculation of C2, as we must calculate X(C, Cp) 
only for C ^ Jl, and not for all C having non-zero fi(C). To calculate the matrix 
element A(C, Cp) for C ^ Jl, we can use the expression (|3"Tj) derived in the previous 
section. All possible C ^ Jl having non-zero /i(C) will either have a particle at site 
fci ~ 1, . . . , L — N — 1 and a hole at site = L, or particles both at k\ = 1 and L with 
a hole at site &2 = L — N + 1, — 1. Let's first consider configurations C £ Jl 

such that the particle outside the queue is placed at fci = 1, . . . , £ — N — 1 leaving a 
hole at site L. If fci ^ L — N — 1 (meaning that there is at least one hole in front of 
the particle), then Aq(C) = 2 for any C" visited in going from C to Cp. The matrix 
element A(C, Cp) for such C is therefore given by 

X(C,C P )=(^\ , h = l,...,L-N-l, k 2 = L. 

If kx = L — N, then Aq(C) = 1 and the rest of the configurations in going from C to 
C P have Aq(C) = 2. This gives 

A(C,Cp)=(-J , h = L-N, k 2 = L. 
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Similarly, in going from configurations with k\ = 1 and k 2 = L — N+l,...,L — Ito 
Cp, Aq is always 2 and therefore X(C, Cp) is given by 



j N L-fc 2 + l 



\{C,Cp)=\^-) , h = l, k 2 =L-N + l,...,L-l. 

If fei = 1 and k 2 = L — N, then Aq(C) = 1 and so A(C, Cp) is given by 



X(C,Cp)=[-) , fci = l, k 2 = L-N. 



Summing all four contributions gives 



C2 



L-JV-1 

E 



fei 



L-N-l 



L — N 



k 2 : 

L-N-l 



L-l 

£ 

=L-Af+1 



L-fe 2 + l 



3 

~~ ~ 2' 

Our method thus gives us ji(r) up to 0(r 3 ) 



j L (r) = r- -r 2 + 0(r 3 ), r < 1, 

which was first calculated by Janowsky and Lebowitz |26j . Unfortunately, this is as 
far as we can go without much effort. To calculate the next-order terms, we would 
have to explore paths starting from configurations having either two particles outside 
the queue and a hole inside the queue, or one particle outside the queue and two holes 
inside the queue. However, tracking movement of two particles or two holes is no 
longer trivial because of the exclusion, and therefore it becomes increasingly difficult, 
albeit possible, to calculate A(C, C). 



Open boundary conditions. Now let's consider open boundaries case with slow site 
placed at site k. Here we can expand the current ji,(a,P,r) in any of the hopping 
rates a, /3 or r. The simplest case to consider is when none of the two remaining 
hopping rates are equal to the one that we are expanding in. In that case B has only 
one configuration: an empty chain if expanding around a = 0, a fully occupied chain 
if expanding around (3 = and a semi-full chain with particles behind the slow site if 
expanding around r = 0. When expanding jh around r — 0, the calculation is similar 
to the one for the periodic case and in fact gives the same result, 

3 

j L (a,(3,r) =r - -r 2 + 0{r 3 ), r<Ca,/3. 

That ji, does not depend on a nor j3 in the small-r limit was recognized long time 
ago by Janowsky and Lebowitz by studying a = /3 case. Note also that in the 
small r limit does not depend on the position k of the slow site either. 
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Around a — or j3 — 0, the calculation is even simpler. There is only one blocked 
configuration with respect to a, and that is an empty chain. This gives us immediately 
a o = ^0 = L i-G- Ci = 1. The expression for c-i reads 

C£J 

where Jq = {C\t\{C) — 1}. There is only one configuration C € Jo that has 
/i(C) ^ 0, and that is the configuration with a particle at site 1, which has /i(C) = 1. 
A similar calculation can be made for the expansion around j3 = 0. For the first two 
coefficient, the final result is thus the same as for the pure TASEP, 



j L (a,P,r) = a - a 2 + 0(a 3 ), a<^,r 
j L (a,P,r) = /3-/3 2 + 0(/3 3 ), /3 « a, r 

Again we see that the coefficients are pure numbers and do not depend on other 
hopping rates, nor on the position of the slow site. 

When a or /3 is equal to r we immediately notice that \B\ (the number of elements 
in B) is greater than 1, and therefore solving (|35|) cannot be avoided. It is the same 
difficulty that we are going to encounter when dealing with more than one slow site 
in the following sections. Let's consider a = r <C (3 case first. Blocked configurations 
can be described as having a queue behind the slow site and an empty segment in 
front of it. Compared to the r <^ a,/3 case, the queue is now no longer of size k (i.e. 
occupying the whole segment behind the slow site), but can be of any size 0, . . . , k 
giving \B\ = k + 1. Let's denote configurations belonging to B with C m , where 
m = 0, . . . , k is the size of the queue behind the slow site placed at k. The equations 
for /o(C m ) are easily generated using (l2Tj) giving 



/o(C ) =/o(Ci) 

2/o(Cm) = /o(C m _i) + /o(C(m + 1)), m = 1, . . . , k - 1 
/o(Cfe) = /o(Cfc-i). 

According to our interpretation using partial exclusion process, this corresponds to 
having a single site with capacity k which exchanges particles with two reservoirs at 
rates a — (3 — 1. The solution to the system above is simply /o(C m ) = const.. If we 
choose g(C) — rr^ [1— Tfe+i], then ag — k and bo = k + 1 so that C\ = ao/bo = k/(k + l). 
A much more involved calculation is required to get the second-order terms. Here we 
state the final result leaving the details of this calculation to |Appendix B| 

, „ N k (k- l)(5fe + 8) , , o, 

j L (a, 0, r) = - 1 ^ 1)2 V + 0(r 3 ), a = r « /?. (42) 

How well truncating the expression (|42[) at second-order approximates j'lM is 
presented in figure 5 where (42) is compared to jl{t) obtained from Monte Carlo 
simulations on a lattice of L = 1000 sites for /3 = 1 and (a) k = 3 and (b) k = 10. 

Using the particle-hole symmetry t, H 1- fL-i+i, a «-> /3, fc <H> L — fc, we can 
also get the expansion for j3 — r <C a which reads 

^^ = L-kTi r -- 4(L- fc + l) 2 ^ 2 + °^)- ^ = »■<«• ( 43 ) 




Figure 5. Current ji^(r,a, /3) as a function of a = r obtained by Monte Carlo 

simulations ( ) on a lattice of L = 1000 sites with /3 = 1, compared to the 

expression i ]4_ > l ) truncated at the second-order ( ) for (a) k = 3 and (b) k = 10. 



In the most complicated case when a = p = r, which corresponds to the partial 
exclusion process with two sites having capacities k and L — k, unfortunately we were 
not able to find even fo(C), i.e. to solve (|2Tj) for general k. This already clearly 
demonstrates that severe difficulties are to be expected whenever \B\ is not small. 

4-2. TASEP with two slow sites 

Periodic boundary conditions. We next consider two slow sites placed at k\ = L — d 
and /c2 = L on a ring of L sites and N particles. We will further assume that 
d < N < L— 43, which means the number of particles in the segment i = L—d+l, . . . , L 
can take values 0, . . . , d. The number of blocked configurations is then \B\ = d+1. The 
corresponding partial exclusion process consists of two sites with periodic boundary 
conditions, which has a simple steady state with all fo{C) = const.. Choosing 
again ji, = (rrfc 1 (C)[l — Tfc 1+ i(C)]) gives immediately ao = d and bo = d+1, i.e. 
ci = d/(d + 1), as in (JT)). The calculation of the second-order term is very similar to 
the a = r case with a single slow site. The final result for jl{t) up to 0(r 3 ) is 

• < d (<f- l)(3d + 4) 2 _ 3 , , tA . 

J L (r) = - 2( ;; 1)2 V + 0(,3). (44) 

Notice that as d — » oo, joo(r) = r — 3r 2 /2 + 0(r 3 ), as in the TASEP with a single slow 
site. A comparison of (|44p. truncated at the second-order, with j'i,(r) obtained from 
Monte Carlo simulations (L — 1000, N — 500) is presented in figure [5] for (a) d = 3 
and (b) d = 10. 



Open boundary conditions. For r <S a, /3, the result is the same as in (|44p . 

■ , d (d- l)(3d + 4) 2 3 , a , Ar , 

jL(r) "d+T r " 2(d+l) 2 r +°( r )> ( 45 ) 

Other cases, i.e. when a or /J are equal to r, are more difficult to deal with as the 
analytical solution to (|2ip is generally not known. 

X Other values of N can be explored as well, but we are here mainly interested in large N and small 
d. 
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Figure 6. Current ji,{r) as a function of r obtained by Monte Carlo simulations 

( ) on a ring of L = 1000 sites and N = 500 particles, compared to the 

expression j44li truncated at the second-order ( ) for (a) d = 3 and (b) d = 10. 



4-3. TASEP with a bottleneck 

Finally, we mention the case of all slow sites being clustered in a bottleneck, which was 
previously studied in [19]. Here the set B is 2^, where I is the number of slow sites. 
Using the previously developed mapping to the partial exclusion process in the limit 
r — > 0, the periodic boundaries case becomes equivalent to the pure TASEP with one 
large but finite reservoir. If we further assume that I < N < L — I so that the finite 
reservoir is never empty nor fully occupied, fo{C) is given by the exact solution of 
the pure TASEP of size I — 1 with open boundaries and a — (5 — 1 . Using the known 
solution of the pure TASEP with open boundaries 06], the current ji(r) up to 0(r 2 ) 
reads 

h{r) = ^±r + 0{r 2 ). (46) 

The same result applies to the open boundaries case for r <C ce, (3, which was 
conjecturecjjj in [T!5]. Unfortunately, due to large \B\ we were not able to find the 
second-order term for general I. (The I — 2 case is already covered by the previous 
example when the distance between two slow sites is 1.) 



5. Further applications 

The approach developed in section[3]is general and can be applied to any unidirectional 
driven diffusive system in which one can identified blocked configurations with respect 
to one of its hopping rates. (Here the unidirectional hopping is necessary to relate 
A(C, C") to weighted backwards paths in the configuration space.) The success of this 
approach will mostly depend on our ability to solve (|2"T]) and ([35]) . If B = {Cp}, these 
are trivially solved and the main problem is to find A(C, C"). For \B\ > 1, we may 
try to solve (|2T]) analytically for small values of \B\ or on a computer for larger values 
using the analogy with the partial exclusion process. 

As a further application that goes beyond binary disorder discussed so far, here 
we mention some results for the slow hopping rates that are not necessary all equal. 

§ In the periodic boundaries case that is true provided I < N, which is the case we consider. 
| Although the notion that the bottleneck behaves as a small TASEP within a big one is not surprising 
and new, it is unclear to us whether the authors of 1191 were actually aware that this picture is exact 
in the limit r — > 0. 
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This type of disorder was in mind in the original idea of MacDonald et al [2][3] , who 
introduced the TASEP to model the process of translation in protein biosynthesis. In 
a simplified description of translation, a ribosome binds to mRNA and moves along it 
codon by codon translating the mRNA sequence into sequence of specific amino acids. 
At each step, the corresponding amino acid is transported to the ribosome by tRNA. 
The availability (abundance) of tRNA is thus believed to be responsible for the time 
scale on which the ribosome moves along the mRNA. Codons with lower concentrations 
of corresponding tRNA will locally suppress ribosome motion across them, acting thus 
as slow sites. An important question, explored extensively in [53], is how are protein 
production rates correlated with specific sequences of codons. Translated into the 
TASEP, the question is to determine the limiting factor for the current with respect 
to the strength and the positions of slow sites. Here we discuss some immediate results 
that stem from our approach applied to the non-binary disorder. We will not include 
another important ingredient for modelling translation - the fact that ribosomes bind 
to approximately 12 codon sites - as it become technically difficult to do it in our 
approach due to the exclusion. 

For two slow sites with hopping rates n ^ r 2 our approach readily gives 

3 

jh{r) = r* - -r* 2 + 0(r* 3 ), r* = min{ri, r 2 } < a, (3 (47) 

This result can be easily generalized to arbitrary number M of slow sites provided 
they all have different hopping rates 

3 

j L (r) = r* - -r* 2 + 0(r* 3 ), r* = min^i = 1, . . . , M} < a, /3 

nj^Tj, i^j 

If two slow sites however share the same hopping rate r, then (j45|) applies provided 
all other slow hopping rates (including a and /?) are mutually different and not equal 
to r. What this tells us generally is that the low-current regime of the TASEP with 
sitewise disorder depends only on the current-minimizing subset of slow sites with equal 
hopping rates, regardless of other slow sites. 

Here we make a modest attempt to test this idea by simulating the TASEP with 
15 randomly distributed slow sites of which 10 have rates r\ (type 1) and 5 have 
rates r 2 > ri (type 2) on a lattice of 1000 sites with a — (3 — 1. Figure [7] compares 
current ji(ri, r 2 , a, (3) as a function of r\ for two values of r 2 , one with r 2 = 0.3 (both 
types present) and the other with r 2 = 1 (only type 1 present). Our data shows no 
significant difference between these two currents for small r\ < r 2 = 0.3. Notice also a 
gap between the current with type 2 slow sites only (- -) and the point r\ = r 2 = 0.3 
where type 1 slow sites turn into type 2, which lowers the current as there are now 15 
slow sites of type 2 instead of 5. 

While this is far from a systematic study, it shows that the low-current regime 
can be safely approximated by the TASEP with binary disorder provided the current- 
minimizing set of slow sites with equal rates is properly identified. It would be 
interesting to test this idea further using real data (e.g. relative abundances of tRNA 
measured for E. Coli 34J) and compare it to the approximate theory of estimating 
currents for the disordered TASEP developed in [55] . 




Figure 7. Current jl(tx, T2, ct, 0) obtained by Monte Carlo simulations on a 
lattice of L = 1000 sites with 15 slow sites of two types, 10 slow sites with rates 
ri (type 1) and 5 slow sites with rates T2 (type 2), plotted as a function of n for 

fixed T2 = 0.3 ( ) and T2 = 1 (O ). Setting T2 = 1 in the latter case means 

that only type 1 slow sites are present. Broken line ( ) corresponds to r% = 1 

and T2 = 0.3 (type 2 slow sites only). In all cases a. = j3 = 1. Slow sites of type 1 
were placed at i = 27, 79, 262, 558, 563, 655, 686, 701, 720, 888 and those of type 2 
at i = 156, 191, 407, 744, 997. 



6. Conclusion 

The matrix-product ansatz and mean-field approximation are both powerful ana- 
lytical approaches for studying driven diffusive systems, the ASEP in particular. 
Unfortunately, in some cases such as site-wise disorder, it is not known how to 
apply the matrix-product ansatz and mean-field approximation is mostly reduced to 
numerical studies. In this article we showed how to access some exact steady-state 
properties of the TASEP with site-wise disorder in the low-current regime, i.e. when 
one of the hopping rates is small. 

Our approach is based on a simple fact, proved in | Appendix A| that the steady- 
state (non-normalized) weights are polynomials in hopping rates. Using this fact the 
steady-state average of any physical observable (e.g. current) can be expanded in 
one of the small hopping rates with coefficients that obey a specific set of equations. 
While our approach is not restricted particularly to the TASEP, it requires that we can 
identify what we call blocked configurations (configurations that the system freezes 
into if one of the hopping rates is set to zero). In that case the zeroth-order coefficients 
are non-zero only for blocked configurations, which drastically reduces the number 
of unknowns. For the TASEP with binary disorder where all slow sites share the 
same hopping rate r < 1, we show that the zeroth-order coefficients in the small 
r expansion are in fact steady-state weights of another (but rarely studied) process 
called partial exclusion process in which more than one particle per site is allowed. 
This mapping, which is exact in the limit r — > 0, is made by replacing all slow sites 
with boxes of capacities equal to the distances between neighbouring slow sites. A 
simple interpretation of this result is that the limit r — > creates a huge separation 
of time scales between particles hopping at rates r and 1, so that any particle that 
jumps across a slow site immediately joins a queue in the front. It is remarkable (and 
discouraging at the same time) that the TASEP with binary disorder even in this 
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simplified case maps to a process which itself is a hard and unsolved problem. 

In cases when zeroth-order terms can still be found, we provide a recipe for 
calculating first-order terms, which comes down to tracking paths in the configuration 
space on the underlying stochastic network. In section U we calculated zeroth-order 
and in some cases first-order terms in the small r expansion of the steady-state 
current for particular disorder configurations previously studied in [19 1 126 ] . As the 
number of blocked configurations increases it becomes increasingly difficult to follow 
this programme analytically. However, because of the mapping to the partial exclusion 
process the reduction of the unknowns is still huge and allows us potentially to use a 
computer instead, even for large lattices. 

Our approach can readily be applied to non-binary disorder with more than one 
type of slow rates, which is relevant for modelling protein synthesis. Remarkably, 
the lowest-order coefficients we get by expanding in one of the slow hopping rates 
do not depend on the other hopping rates. In other words, the low-current regime of 
the TASEP with site- wise disorder depends only on the current-minimizing set of slow 
sites with equal hopping rates, regardless of other slow sites. Once this subset is found 
(which may be a hard problem in itself), one can work with binary disorder only and 
go from there using either the approach developed here or using e.g. phenomenological 
approach that looks for the largest cluster of slow sites [27j[28] . 

Finally, we mention that our approach can be applied to any other driven diffusive 
system with particles hopping unidirectionally, provided we can identify blocked 
configurations with respect to one of the model's hopping rates. In the network 
language, unidirectional hopping is necessary to avoid loops that may exist even if 
one of the hopping rates is equal to zero. For example, our approach does not work 
for Langmuir kinetics |20j or a multi-lane exclusion process where particles can change 
lanes [22] . However, it could be useful for studying more complex lattice geometries 
without resorting to mean-field approximation or even to understand why in some 
instances the mean-field approximation is satisfactory. 
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Appendix A. A formal solution of the stationary master equation 

Consider an ergodic continuous-time Markov jump process with transition rates 
W(i —tj)- Master equation for the steady-state probabilities Pi is then given by 



where M is the total number of states. Equation (|A.1[) can be written in a more 
compact form by introducing a matrix L^, 




(A.1) 




(A.2) 



j 




(A.3) 



k 
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Note that L is a left stochastic matrix, which means that Y]. Ljj = for any j. 
An important property of a left stochastic matrix is that it has all of its cofactors 
Cij = (— 1) 4J Mij independent of i. (Here My is defined to be the determinant 
of the submatrix of L obtained by removing i-th row and j-th column from L.) 
Since we assumed that the process is ergodic, the stationary state satisfying equation 
y\. LijPj = is unique and non-trivial, which means that detL = 0. By using the 
Laplace expansion for detL and the aforementioned fact that Cjk = Ckk we get, 

det L = = E L ikCjk E L ^C kk , (A.4) 
fc k 

which means that Pi must be proportional to Cu. Since C'u is the determinant of 
a matrix whose matrix elements are linear combinations of the transition rates, we 
conclude that Pi must be a polynomial in all the transitions rates present in (IA.1I) . 



Appendix B. Second-order term C2 for a = r in the TASEP with a single 
slow site 



The starting point is the expression for ci 



ai bi ci ci - 1 

c-2 = T ci— = (ai - oij— + ai — 

bo b b b 

= -<^S/ l(c,+ <^c5/ i(C) ' 

where Jk — {C\rk(C) — l,Tk+i{C) = 0}. Summations in (|B.1[) can be further 
separated into 



CiJ k c« t c{j k 

C£B C£B 



E f^ G )= E E w)- ( B - 2 ) 



C£j k cej k cej k 

C€B C&B 



The first sum in (jB.lj) can be greatly simplified by that fact that there is only one 
C G B for which g(C) = rrk(C)[l — Tk+i(C)] = 0, and that is an empty lattice, 

E /i(C) = /i(Co)- 

C<tJ k 
CEB 

To find /i(C) for C £ B, we have to solve ([55)1 , which in this case reads 
h(C ) =/i(C 1 ) + / ll (Co) 

2/i(C m )=/ 1 (C m _ 1 )+/ 1 (C m+1 ) + /i 1 (C , m ), m = l,...,fc-l (B.3) 

h(c k ) =/i(c fc _i) + fti(c fc ), 

where hi(C m ) is given by (|36|) . Luckily, this system admits closed expression for 
h(C m ) which is 

fi(C m ) = /i(C ) - E* ' hi{C m -i), m = l,...,k. 
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The first sum in (IB. 21) can be therefore written as 



E /i(c) = fc-/i(Co)-ED-MC m -i) 



ceB 



m—1 i—1 

2 



fc-l 



i=0 



Altogether, the expression for C2 can be written as 



C-2 



(fc + l) 2 



fc-l 



C$B 



E 

C(fB 



(B.4) 



Now, let's go back to the expression (|36|) for h±(C). For a given C m , the set 
Si(C m ) n £? consists of all configurations C having /i(C) ^ that can be reached 
from C m by moving particles backwards and hitting a slow edge in the last move only. 
Starting from C m (with m particles in the segment i = 1, . . . , k by the definition), the 
resulting C will have either m — 1 or m + 1 particles in the same segment depending 
on whether we crossed the left boundary or site k in the last move, respectively. The 
resulting C will necessarily have a hole at site 1 in the former case and a particle-hole 
par at k, k + 1 in the latter case. Thus by defining J\ = {C\t\{C) = 0} in the same 
spirit as we defined Jk, we can write 



E f^= E h(c)-i{J2^(c) = i} 

iSi(c ) cej k i—\ 

C(B C£B 

k 

E = E • ME^cq = m - i} 

:Si(c m ) cej i=l 

fc 

+ ^/ 1 (C)-l{^T,(C)= m + l}, m = l,.. 

cej t »=1 

fc 

E E fi(c)-i{Y,n(c) = k-2} 



cesi(c fc _ 1 ) cej ,:=i 

C£B C£B 

fc 

E = E h{C) ■ l{J2n(C) = k 1} 

cesi(c fc ) cejo j=i 

C(B C(B 

where 1{X = x} is an indicator function, 1{X = x} = if X ^ and 1 if X = x. 

Similarly, V\{C m ) in (|36|) consists of all configurations C having /i(C) 7^ and 
precisely m particles in the segment i = 1, . . . , k. We can further divide this set into 
subsets depending on whether C belongs to Jo D Jk (for which A r {C) — 2), J but 
not Jk and vice versa (for which A r (C) = 1) or to none of these sets (for which 
A r {C) — 0). The second sum in (1551) can be then written as 
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k 



E m°) = E /i(c , )-i{E r ^ c ') = } 

y l( c ) «g 

k 

E /i(c) = E f^°) ■ ^E^) = w 

E /i(C)-l£r,(C) =m}, m=l,...,fc-l. 



C£B 



To ease the notation let's introduce a^\m) and a[ k \m) defined as 



a ( °\m) = E /i(C)-l£r,(C) = m}, m = 0, . . . , fc, (B.5) 



cej i=l 

C£B 



= E • 1{E T< ( C ) = m >' "» = 0,...,fc. (B.6) 

ceJt i=l 
a$B 

From here it is easy to see that 

4 0) (fc)=0, 4 fc) (0) =4 fc) (fc) =0. 
We can now write h\(C m ) as 
/H(C ) =4 fe) (l)-4 0) (0) 

hx(C m ) =a ( f ) (m + l)-a i f ) {m) + a i °\m-l)-a i ° ) (m), m = l,...,fc-2, (B.7) 

= a< 0) (fc - 2) - a' 0) (fc - 1) - a[ k \k - 1) 
Using definition for a[ k \m) we can also write 

k-l 

E A(c) = E a i fe V)- (B-8) 

Cej k m =l 
cgB 

While each af (m) and (m) can be calculated explicitly using (l3~Tj) . we can 
further simplify calculation by showing that we must only calculate their difference 
a^\m) — a[ k \m). Using (|B.7[) . we can show after some algebra that 

1 fc-i fe-i 

2 E( fc - *)( fc - * + ^mco = -fc«i 0) (o) + E a i fc) ( m ) 

+ (fc — m) [a [ k ^ (m) — (m)] . 

m— 1 

Using (f3~Tj) it is also straightforward to calculate the second sum in C2 . Depending on 
the number of particles m in the segment i — 1, . . . , k we can show that 

E /i(<7)=a< 0) (0) + |(A-l). (B.10) 

C^B 
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Inserting (|R9|) . (jR8|) and (|B~TOl) in (|R4|) gives finally 




fc-1 



k(k — 1) + (fc — w) ai (to) — (to) 



fc > 2. 



The difference aj ' (to) — (to) can be further written as 

k 

a?\m)-a<* ) (m) = £ ji(C) ■ 1{C| 5>(C) = m} 



i=l 



(B.H) 



^ / 1 (C)-l{C|^r i (C') = m}. 



C(B 



Using (f5T|) we can easily find that 



fc-i 

E E 

m=l cej ,cia k 

G4B 



(fc 



fc 

m)/i(C).l{C|5>(<7) 

i=l 



TO,} 



(fc-2)(7fc- 1) 



fc-i 



To compute the second sum in (|B.III) we must first identify configurations C ^ J7q 
and C £ Jk that give /i(C) 0. Since sites 1 and fc must be occupied and site fc + I 
must be empty, we have only one option for m = 2, . . . , fc — 2 and two options for 
to = fc — 1. For m — 2, . . . ,k — 2, the only way we can reach a blocked configuration 
by moving particle backwards is to start from a configuration that has a queue behind 
the slow site and a particle at site I. The blocked configuration with m — 1 particles 
is then reached by moving a particle at site 1 backwards, which gives A(C, C m -i) = f . 
We can thus write 



fc-i 



fc-i 



J2 E fi(C)-l{C\J2n(C) = m}=J2(k-m)-l 

i=X rn—2 
k 

+ h(C)-l{C\Y / r l (C) = k-l} 



m=l C<tJ ,C&J k 
C<£B 



(B.12) 



C<tJ ,C£J k 



i=l 



For to = fc — I there is a single hole in the otherwise full segment i = 1, . . . , fc, so in 
addition to moving backwards a particle at site 1 we can reach a blocked configuration 
with to = fc particles by moving a particle from the segment i = fc + I,...,Lor from 
the right reservoir. Computing A(C, C&) is then straightforward but more complicated, 
because we must count all the ways in which a particle and a hole can both move until 
Cfc is reached. Instead of calculating each A(C, Ck) individually and then make the 
summation, we will use the following trick. Let's denote with C^.^ configurations 
that have a hole placed at fci = 2, . . . , fc— 1 and a particle placed at ki = 2, . . . , L—k+1, 
where both k\ and fc2 are measured relative to the site fc. Here fc2 = L — fc+ I denotes 
a particle in the right reservoir, i.e. a configuration with no particles in the segment 
i = k + 1, . . . , L. For a given ki = 2, . . . , k — 1 the equations for j \(Gk\,k^) read 
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2/i(Cfc 1 ,fc 2 ) = fi(C kl -i,k 2 ) + h{Ck x ,k 2 -i), k 2 = 2,...,L-k-l 

(1 + P)A(C kl , L -k) = fi{C kl - 1M ) + h{C kl , L -k-i) 

fl(Ck u L-k+l) = h{Ck!-l,L-k+l) + PfliCkx^-k) 

Summing all the equations for a fixed k\ = 2 . . . , k — 1 gives after some algebra 

L-k+l L-k+1 

E /i( c 'fci,fa)= E /i(c fcl -i >fe2 )+/i(fci,i)- 

fc 2 =2 fc 2 =2 

The advantage here is that we is no dependence on j3. Solving this is as a recursion 
relation in fci we get 

L-fc+l L-k+1 fci 

E h( C MM)= E /i(Ci,fc 2 )+E/i( C M) 

fc 2 =2 fc 2 =2 i=2 

Both sums on the r.h.s. are now simple to calculate using (|3ip because either particle 
or hole is always fixed. The final result is 

fc-l L-k+l , v fc-1 

E E mx,fa) = *-2 + u ■ (R13) 

fci fc 2 =2 ^ ' 

Inserting (|B.13|) in (|B.12I) and then in the expression for C2 we finally get 

(k- l)(5jfc + 8) 
C2 = - 4( fc + l) 2 ' ^ L 
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